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Sources of isocurvature perturbations ar\d large non-Gaussianities include field degrees of free- 
dom whose vacuum expectation values are smaller than the expansion rate of inflation. The inho- 
mogeneities in the energy density of such fields are quadratic in the fields to leading order in the 
inhomogeneity expansion. Although it is often assumed that such isocurvature perturbations and 
inflaton-driven curvature perturbations are uncorrelated, this is not obvious from a direct computa- 
tional point of view due to the form of the minimal gravitational interactions. We thus compute the 
irreducible gravitational contributions to the quadratic isocurvature-curvature cross-correlation. We 
find a small but non-decaying cross-correlation, which in principle serves as a measurable prediction 
of this large class of isocurvature perturbations. We apply our cross-correlation result to two dark 
matter isocurvature perturbation scenarios: QCD axions and WIMPZILLAs. On the technical side, 
we utilize a gravitational Ward identity in a novel manner to demonstrate the gauge invariance of the 
computation. Furthermore, the detailed computation is interpreted in terms of a soft-^ theorem and 
a gravitational Ward identity. Finally, we also identify explicitly all the counterterms that are neces- 
sary for renormalizing the isocurvature perturbation composite operator in inflationary cosmological 
backgrounds. 
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I. INTRODUCTION 



As physics beyond the Standard Model is expected to contain many fields in addition to the inflaton, there are 
many candidates for isocurvature perturbations in the context of inflationary cosmology including those of the dark 
matter Indeed, the current data is consistent with the existence of an 0(5%) isocurvature component [1-8|. Fur- 
thermore, it is well known that quadratic isocurvature perturbations (i.e. the vacuum expectation value of the field 
is much smaller than the Hubble expansion rate) are one of the very few ways to generate measurably large local 
non-Gaussianities |9 -40| in the context of the slow-roll inflationary paradigm. The only nontrivial requirement that 
the isocurvature field degree of freedom must possess is that it be light enough to be excited by the inflationary 
quasi-de Sitter (dS) background and that it not be conformally invariant. In the literature f4r-'43l, quadratic isocur- 
vature perturbations are often assumed to have negligible cross-correlations with the curvature perturbations (which 
corresponds to the inflaton field degree of freedom dressed by gravity). However, the gravitational interactions lead 
to a minimum cross-correlation, which in principle can be observationally important. We present a computation of 
this minimal gravitational cross-correlation in this paper. 

As explained below, the form of the gravitational interaction between the curvature and isocurvature perturba- 
tions naively suggests that there can be cross correlators which do not vanish in the long wavelength limit. If this 
was true, the cross correlation can dominate over the isocurvature two-point function in the observables since the 
latter vanishes in the long wavelength limit for a massive field. By an explicit rigorous computation, we show that 
the cross correlator vanishes in the long wavelength in such a way that the cross correlation induced by gravity never 
dominates over the isocurvature two-pornt function, given that the curvature inhomogeneity perturbation is char- 
acterized by a strength of order 10^^. We explain this qualitatively as well using a combination of a soft-^ theorem 
I44--61 j and a Ward identity associated with a spatial dilatation diffeomorphism. We also check the gauge invariance 
of our computation using a Ward identity. 

Among the possible isocurvature candidates, thermal dark matter is usually produced copiously by the inflaton 
decay products, which typically leads to a large suppression of isocurvature effects. On the other hand, nonthermal 
dark matter that is not produced by the inflaton decay can easily generate large isocurvature effects that survive 
until today. Hence, as an illustration, we apply our computation of the cross correlation to two different nonthermal 
dark matter models: QCD axions and WIMPZILLAs. In both cases, we find a cross-correlation characterized by 
the parameter |/3| ^ 0(10^^) (the parameter definition is given in Eq. |29 ) which is below the boundary value of 
0(10 "2) when the cross correlation becomes competitive with the isocurvature two-point function. In principle, /3 
can be measured and is a generic prediction of this class of nonthermal dark matter quadratic isocurvature models. 
Note that even though the nonthermal dark matter fields can be identified with the isocurvature degrees of freedom, 
this scenario is consistent with the WIMP dark matter scenario since the isocurvature perturbations can be as small 
as an order 10^^ fraction of the total dark matter and still leave an isocurvature imprint on the CMB spectrum. 

The order of presentation is as follows. In Section|llJ we present our assumptions about the inflationary cosmology, 
review gauge invariant variables in the perturbation theory, and summarize the observational constraints on the 
isocurvature scenario relevant to our paper. One of the most important aspects of this section is our review of 
features of the j6 variable that we compute. In Section III we first explain two naive estimates, one leading to the 
wrong observationally large result, and the other leading to the correct suppressed result. In explaining the correct 
estimate (which requires assumptions that cannot be known without the justification of a full computation), we 
present the interpretation in terms of a soft-^ theorem and a Ward identity. The rigorous explicit computation at 
one loop is then presented, demonstrating how the correct naive estimate result is achieved. We also present in this 
section how gauge invariance is achieved for these quadratic isocuvature computations using a gravitational Ward 
identity. Next, we apply these results to the axion and the WIMPZILLA scenarios in Section IV This section contains 
a detailed explanation for choosing nonthermal dark matter to illustrate the computations of our paper instead of 
thermal dark matter. Finally, we summarize our results in Section |V] In appendices, we collect technical details and 
also supplementary computational results: the radiation transfer functions is derived in Appendix [A] a brief review 
of the gravitational Ward identity used for the gauge invariance computation is given in AppeiiHix |B| the ADM 
formalism is reviewed in Appendix |C| the details about the Pauli-Villars regulator is explained in Appendix [P] and 
the two point function computation in the uniform curvature gauge is presented in Appendix |E] 



II. A CLASS OF CURVATURE AND ISOCURVATURE PERTURBATIONS 



Inflation through quantum correlator djmamics generates "classical" initial conditions for superhorizon cosmolog- 
ical fluid perturbations |62-65|. The resulting initial conditions for the classical equations governing classical fluid 
variables (which are set during radiation domination before the CMB last scattering time) are categorized into two 
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t3^es: adiabatic and isocurvature f66'-'69l. An adiabatic initial condition is intuitively characterized by all species 
composing the fluid having the same initial number overdensities. In the context of inflation, if there is a single dy- 
namical degree of freedom (p during inflation such that after a few efolds of inflation, the quantum vacuum boundary 
can be approximated as Bunch-Davies initial conditions (for a discussion of number of efold requirement see e.g. 
f70 1), and if all the degrees of freedom during radiation domination come from the inflaton decay, then this adiabatic 
condition is the resulting approximate classical boundary condition during radiation domination era of the universe. 
An isocurvature initial condition intuitively corresponds to setting nonzero the initial difference of the number over- 
densities of at least one pair of fluid element species while setting to zero the total energy density inhomogeneity 
on long wavelength scales. Because these two types of initial conditions are linearly independent, a generic initial 
condition to the linearized perturbation equations can be written as a linear combination of them. 

In this paper, we are concerned with the following physical system which is generic for isocurvature scenarios. 
One real scalar slow-roll inflaton degree of freedom (p dominates the energy density during inflation. During this 
time period, there exists also another light degree of freedom a which has no coupling to (p stronger than gravity. 
We assume that this system carries an approximately conserved discrete charge (such as Z2 broken at most by 
a model dependent non-renormalizable operator) such that the one particle states are stable and can act as dark 
matter. Note that since we do not require all of the dark matter to come from a, this system is consistent with the 
existence of the weakly interacting massive particle (WIMP) dark matter. If WIMP dark matter exists, the parameter 
cOcr = Cicr /CicDM < 1 will play a role, and this scenario can yield interesting isocurvature signatures for cOcr as small 
as 10^^ 1.30.1 . The action of this system can thus be written as 

S[cp.cr,W]= I (dx) I ^M^R + [- ^g^^'d^cpdy^p -V{cp)] + [- \g^^^^a^ya - !J(£r)] | + S,,, \<p, {^p}\ (1) 



where K is the Ricci scalar, Mp = (dx) = d'^x^\ det(^;,i,)|/ and S^;, corresponds to the action of the reheating 

degrees of freedom {)/'}. We assume that {i/'} is heavy during inflation such that it can be integrated out or if {i/;} are 
light, they are conformal such that they are not excited during inflation. After inflation ends, we assume {ip} fields 
are light, leading to a successful reheating scenario. The only special initial condition dependent assumption that we 
make in this isocurvature scenario is that {a) <^ H / {2tz) during inflation even when d^U{(7)/da''^ <^ H. ^ Because 
{</) = during inflation, a by itself does not spontaneously break time translation invariance and therefore does 
not mix with S(p in forming the gauged time translation Nambu-Goldstone boson Hence, we can treat the scalar 
fluid variable [{Sg^y, 6(p) as the curvature degree of freedom and [a, Q as the isocurvature degree of freedom. (As 
we will show in detail below, the isocurvature degree of freedom will be quadratic in a and will involve ^ as a 
difference). 

Thus, the basic physics picture of the classical fluid that we are concerned with in this paper is the following. 
To predict CMB temperature fluctuation (ATAT), we must compute the cross correlation (^sC) since at the lin- 
earized level, Einstein-Boltzmann equations give the relationship AT /T ^ Ci^ + €2^$ for computable order unity 
(for long wavelengths) coefficients c,. Up until this paper, there has never been an explicit computation of the 
{SsQ I \/ {^0 i^s^s) coming from irreducible gravitational interactions.^ What will emerge is a clean universal result 
that applies to a wide range of isocurvature models including those of the QCD axions (in a particular initial condi- 
tion regime) and WIMPZILLAs. We find that (^sC) contribution is generically subdominant to {Sg^s) in the case of 
pure gravitational interactions. 

In the following, we establish our conventions in describing this isocurvature degree of freedom carrying the 
non-adiabatic initial condition information. In the process, we review the gauge invariant construction of these 
cosmological perturbations and the current CMB observational constraint, which represents the strongest constraint 
on the isocurvature initial condition derived from inflation. 



^ Note that even with a Gaussian distributed values of (17) on an inflationary patch with a Gaussian width H/ (27r), there is abouta2/3probability 
that such initial condition configurations can be found. Also, an unbroken discrete symmetry such as Z2 : cr ^ —a can stabilize the VEV. In the 
context of supergravity, generic terms in the effective potential however can appear leading to (u) 7^ during inflation. In the end, whether or 
not (c) = is model dependent, but it is not fine tuned. 

^ As we will later explain, we do not compute {SgQ analytically fully beyond the time of the end of inflation. However, the importance of the 
isocurvature cross correlation can be generically predicted by {SsCl / \/ {^Cl i^s^s) which is insensitive to the post-inflationary evolution for 
superhorizon modes. 
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A. Gauge Invariant Construction 



The cosmological inhomogeneity perturbation variables are generally spacetime coordinate gauge-dependent be- 
cause of the coordinate dependent definition of fictitious background metric slices. From the perspective of matching 
classical equation initial conditions to inflationary quantum correlator computations, identifying gauge invariant 
combinations is helpful |71-73|. On the other hand, the gauge freedom involved in computing gauge invariant 
quantities facilitates the quantum computation. Hence, understanding the gauge dependences of the correlation 
computations is helpful. In this subsection, we review the gauge invariant variable construction and establish our 
notation. For a more general discussion, see for example ll7Tll72l [74-84l. 

is) 

In (f, x) coordinates, we parameterize the metric as = gjtv + ^gfiv where the scalar metric perturbation is 



the background metric is g^y = diag{ — l,a-^{t),a^{t),a^{t)}, and derivatives are denoted as usual as X , = dX/dx'. 
Under the diffeomorphism x ^ x + e where 

ef = ie',a'me')), (3) 
the scalar metric perturbation components transform as 

AA = -2He°, AB = -4e^ (4) 
A£ = -2e°, AF =^{e° -e^ +2He^) (5) 

which is obtained from Sg^fJ SgfJ + AiSgfJ) with AiSgfJ) = - C,e)'d^,g¥y 
Similarly, we parameterize the perfect fluid stress tensor for a fluid element a as 

tW=tW+Jt(:) (6) 

where T^y = dmg{pa, Pa, Pa, Pa} contains the average energy density and pressure seen by a comoving observer, 
ST^.^ = PaSgf + a^&ij&Pa, Srjo^ = P(a)Sgio^ " ipa + Pa)Sul"^ (where SUI"^ is the velocity perturbation), and ST^l> = 
~Pci^S(m ~^ ^P"' Under the diffeomorphism of Eq. jsj, the energy density perturbation transforms as 

ASpa = -e%. (7) 

In practice, gauge-invariant variables are constructed by combining metric perturbations and other perturbations, 
such as densities. A popular choice is 

^«-y-Hfl. (8) 

2 Pa 

For example, the first-order gauge-invariant perturbation associated with the inflaton (p is usually defined as 



(see for example Ref. Il73l and references therein). Now, one can form a quantity that is conserved through reheating 
by defining 

^tOt = L^rCr (10) 
i 

where 

L^n Pn ~v ^ n 
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Because there must be reheating d3mamical degrees of freedom, ^^q^ must involve at least 2 degrees of freedom by 
the end of inflation of any single field slow-roll model. In single field slow-roll scenarios, what is done in practice is 
to argue that the reheating degrees of freedom are integrated out during inflation and then integrated back in at the 
end of inflation due to the different location of the inflaton VEV at the end of inflation. Alternatively, another often 
used assumption is that the main reheating degree of freedom are conformal such that no isocurvature fluctuations 
are appreciably excited during inflation. This means that in single field models, we have 

^tot ~ (12) 

up to ambiguities in how one hides the reheating degrees of freedom. 

One reason why the combination of Eq. | |10| | is convenient is because the superhorizon mode of this is approx- 
imately conserved through reheating if this mode object can be shown to obtain an initial conditions of what is 
sometimes referred to as the adiabatic solution |45, 73 1 and there are no non-adiabatic processes that mix super- 
horizon modes of isocurvature degrees of freedom with ^IqI- Such classical adiabatic solution initial conditions are 
generated by the Bunch-Davies quantum fluctuations for and we will restrict the couplings of the isocurvature 
degrees of freedom (discussed below) such as to avoid non-adiabatic mixing. This means that Eq. | [T2| ensures that 
^tot is approximately conserved if + Pip dominates over others. More explicitly, as discussed in the introduc- 
tion to this section, suppose there exists only one isocurvature field degree of freedom which we call a during the 
inflationary period.^ The total curvature perturbation can be written as 

^tot = ^<P+rcri^^-^'p) (13) 

with the sum over n runs over cp and a (assuming that ip has been integrated out during inflation). However, one 
can estimate that the coefficient of during inflation is 

r. < j^.Al - 10-" (14) 



which makes the approximation of ^^-gt ~ (,(p accurate, just as in the single field case of Eq. (12 1. Thus just as in 
the single field scenarios without a, ^^-q^- acquires an approximately adiabatic boundary condition from the Bimch- 
Davies vacuum field fluctuations. 

To complete the examination of how ^^-q^- is used in the scenario of concern in this paper, let's look at the time 
period surrounding the reheating transition when the universe reaches radiation domination. Near the time of the 
completion of the reheating, the variable ^^ot approximately 

ftot ~''<pC<?) + E''>f;^V/ (15) 

i 

such that after the inflaton decays, we have = and 

^tot~E''ft^*r (16) 



^ (The approximation used in Eq. 1 15 1 neglects the r^r contribution because of Eq. 1 14 .) It is also a standard assump- 
tion that 

= ^tot. (17) 

which is rigorously true if one relativistic species dominate the fluid (e.g. r^j ~ 1) or if the decay process does 
not redistribute the spatial inhomogeneities of ipi in a distinct configuration from that of (p. ^ This justifies the usual 
statement in the literature that ^^gt defined in Eq. | [T0) is primarily useful for arguing how a combination of quantities 
involving the inflaton and the reheating decay products remain unchanged through the reheating phase transition. 
Here, we have merely described how this argument is not changed by the presence of a because of the smallness of 



Tg- in Eq. ( 14 1 during the primordial periods of interest. 



^ The species a will later be identified dark matter candidates such as the axions and WlMPZILLAs. 

* In the case that Jpj is integrated back in at the end of inflation, we have made the assumption that this does not change ^jpt 
^ However this need not be true for more general reheating scenarios. 
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In summary, as long as boundary conditions for the classical fluid equation are evaluated at a time when ra is 
small (compared to the accuracy desired), we can neglect the Yg- contribution from ^^-q^- both through reheating and 
until the time that boundary conditions for the classical fluid equations are imposed. Hence, if ^^ot remains constant 
on long wavelengths (due to the initial conditions set by the Bunch-Davies vacuum), Eqs . | [T3| and ( [T4| imply that 
the effective curvature perturbation during this early primordial epoch is given by Eq. < [12| . Hence, in me discussion 
below, we will drop the (p subscript and write 

C = « ^tot- (18) 

During this radiation dominated early primordial time tp, the relationship between super horizon A{tp,k) and the 
value of f (fg. A:) evaluated at the end of inflation time tf, is 

A(t„,k) 2 , 

« (19) 

in the Newtonian gauge (B = F = 0) and the presence of gives a small error controlled by Xa- 

At the same radiation dominated era^ when initial condition is set by ^^-q^- k, ^, the inhomogeneity of the small 
mixture of dark matter component a can be related to the isocurvature perturbation Conventionally, this infor- 
mation is parameterized by the gauge-invariant isocurvature perturbation i67ll69]|76l 

^s(f,fc)^3(^.(f,fc)-^tot(f.^))- (20) 

The physical interpretation of this quantity can be see by noting that when a particles are dominantly non-relativistic 
and the universe is radiation dominated, this expression becomes 

.s(a) = ^-^^ (21) 

where represents the photon energy densities. This clearly represents the difference in number densities of cr 
and '^l Assuming that the radiation inhomogeneity is characterized by ^ as explained in Eqs. (17i and (18 i during 
radiation domination, we have 

^s(f,fc)«3(^^(f,fc)-af.,^)) (22) 

Similarly to the case of ^tot' lo^§ wavelength limit of generated from Bunch-Davies initial conditions simplify 
(partly because of causality) in the absence of non-adiabatic processes mixing of with other superhorizon degrees 
of freedom. The l,a mode for a comoving wave vector k becomes constant once \kl a\ ^ H and ^ H because the 
mode functions involved in Z^^j are governed by the Hubble friction once these conditions are satisfied. 

Although the key correlator computation result of this paper involving /5 evaluated at the end of inflation is in- 
dependent of the transfer function evolving the isocurvature degrees of freedom after the end of inflation, because 
its immediate phenomenological application to CMB requires a transfer function describing this post-inflationary 
evolution, we will restrict our illustration in Section |IV] to the situation when the chemical reaction rates that mix 
a and the radiation co mpon ents are negligible. We will discuss in more detail the cross section constraint for this 



condition in Appendix |IV A 



B. Observational Constraints on Isocurvature Perturbation 



The current observational data shows that the CMB power spectrum is consistent with the adiabatic initial condi- 
tions. However, it does not rule out mixed boundary condition contributions from CDM isocurvature perturbations. 
Schematically, the temperature fluctuations depend Linearly on ^ and 5e, initial conditions as 

^=ciC + C2<5s (23) 



* During this time period, there is possibly a population of thermal dark matter components such as thermal WIMPs. 

^ It is interesting to note that since number densities can diverge while gravitational physics does not care about number densities (in favor of 
energy densities), this choice of variables is unfortunate in situations when there are IR divergences. In this paper, we stick to this convention 
which is prevalent in literature. 
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where c, ^ 0(1). Hence, the CMB temperature correlation data constrains 



fc3 f d^p ,AT(p)Ar(fc), ^^2^^^ 



2712 7 (27r)3^ T T 
where (SI 



(24) 



^(ap)r(^)> = A!w^ (25) 



'^'^(<5s(pMK^))=Alw¥ (26) 



(271)3^^^'^^ '>SW fc3 



--3(<5s(p)rW>=A2,,^(fc)^ (27) 



^ - Kf(^)tip)' ^^^^ 



(29) 



^which are customarily evaluated in the primordial epoch when k corresponds to a far superhorizon scale such that 
the A^(A:) objects are constant in time. Typically the data constraints are parameterized by evaluating a and /5 at a 
pivot scale A: = fcg f3^ 85 1. An important utility of this parameterization is the following fact: a necessary and sufficient 
condition for the cross correlation to be a significant part of the isocurvature contribution is to have |/3| > |c2/ci | 
for a < 1. For example, in order to have approximately the same level of the angular power spectra from both pure 
isocurvature correlation and and cross-correlation at the intermediate scale I ~ 200, i.e. qP"^^^^° ^ Qcrosscor^ 
fractional cross-correlation should satisfy ^ 4 x lO^^ Another utility of the /5 variable comes from the fact that 
when there are non-trivial transfer functions governing A^^^ and A|^ after the end of inflation, the transfer function 
factors can cancel in the expression for /5. We will use this feature later to compute /5 based on just the (quasi)-dS 
mode function behavior^ 

As far as the experimental numbers are concerned, the isocurvature contribution to the CMB temperature per- 
turbation is expected to be roughly less than 10% compared to the curvature contribution. More precisely, the 
Planck+WP limits ^-'S] are 

a|^^0 < 0.016 (95% CL) and < 0.0011 (95% CL), (30) 

where the isocurvature power spectrum is assumed to be scale-invariant, i.e. rijso = 1. The significant difference in 
the upper-bound of a between uncorrelated and totally (anti-)correlated cases can be explained by the ratio j6 / y/a 
already discussed above. The difficulty in improving the current isocurvature bound with data on short wavelengths 
can be seen in Fig. nj where one sees a fall-off of the isocurvature spectrum on short scales (? > 100). This fall-off 



is generic and can be attributed to the transfer function effect encoded by Cj [k)/ C2{k) in Eq. 1 24 1 for k > kgq (where 
keq / ^ 10^2 Mpc~^ is the wave vector associated with matter radiation equality). To understand why Ci (k) / (k) 
generically becomes large for k > ke^, note that isocurvature modes with k > keq enter the horizon during radiation 
domination. Because the isocurvature effect on the temperature spectrum is gravitational, the value of cj (k) / C2{k) is 



Our sign conventions are such that negative values for fi correspond to a positive contribution of the cross-correlation term to the Sachs-Wolfe 
component of the total temperature spectrum. See, e.g., f3'.'851. 
' We will use the exact dS approximation for the massive a and use the quasi-dS approximation for only the massless scenario. The corrections 
coming from the the deviations away from the exact dS background in principle can be absorbed into the transfer function multiplying the 
superhorizon mode function which cancel out in fi due to a common appearance in the mmierator and the denominator. 
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Figure 1: Angular power spectra Q from pure adiabatic(dotted), pure isocurvature(dashed) contributions. The solid red line 
corresponds to be qualitatively the contribution of the absolute magnitude of the third term in Eq. (j24|. The plotted pure adiabatic 
perturbation has the spectral index rig = 0.96. For isocurvature perturbations, the spectral index is 1 and the isocurvature 
fraction a = 0.067 defined at fcg = 0.002 Mpc^^, and the fractional cross-correlation |/3| is 1. 



proportional to the ratio p^{t{k)) / p„-{t{k)) of the radiation energy density to the energy density in the isocurvature 
degree of freedom at the time t{k) when mode k > k^q enters the horizon. Since shorter wavelengths enter the 
horizon earlier, pji{t{k))/ pcr{t{k)) is larger for shorter wavelengths, making (k) / C2{k) larger. For those readers not 
familiar with this physics, some of the details of the transfer function are reviewed in Appendix [A| 

Because of the large differences in the constraints between /5 = and /5 = — 1, estimating trie cross-correlation 
is crucial to restrict parameters and give observable predictions of isocurvature models. In particular, the axion 
scenario with a negligible homogeneous vacuum misalignment angle (and similarly the WIMPZILLA scenario with 
a negligible homogenous backgroimd field value) predicts detectable non-Gaussianity Il24ll30ll86l 

provided the assumption the cross-correlation is zero, i.e. /5 = 0. However, as we will explain, this assumption is 
not obvious for massive field quadratic isocurvature scenarios, and the reexamination of this assumption is one of 
the goals of this paper. 



III. COMPUTATION OF CORRELATORS 



In order to provide the initial condition of the classical fluid equations, it is standard to compute the quantum 
equal time correlators with the inflationary background approximated as a Bunch-Davies vacuum. In this section, 
we compute the correlators using the "in-in" formalism (e.g. see Weinberg [87]). More specifically, in the context of 
canonical quantization, we perturbatively compute the expectation value of an operator Q{t) 



(QW)=E(-o" r dh r dt2- 

n 



dt„ 



(32) 



where the superscript 1 stands for the interaction picture and Q{t) represents a product of canonically quantized 
operators. 

In the scenario explained in Sec. Illj we consider the gravitational coupling whose interaction Hamiltonian is de- 
rived from the ADM formalism witha given choice of gauge. For the computation of the cross-correlation to leading 
order in gravitational coupling, we need at least up to the cubic coupling H^^^tj' where cr is a spectator field during 
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inflation. The interaction Hamiltonian is difleomorphism gauge-dependent. For two commonly used gauges, the 
comovrng gauge{Scp = 0) and the uniform curvature gauge(A = 0), we have 



d'xa'{t)Trit,x)5g^,y{t,x), 



-2i 



(33) 
(34) 

(35) 



where is the stress energy tensor of the field a, and Sg^y is the metric perturbation and the superscript (C) and 
(Lf) denote the comovrng gauge and imiform curvature gauge, respectively. A detailed derivation of the interaction 
Hamiltonian using the ADM formalism is presented in Section [C| 

The isocurvature perturbation So, should be also written in terms of quantum operators associated with the energy 
density pa of the particle a. Since the energy density pa is written in bilinear form of a and since the energy density of 
CDM are often those of non-relativistic particles at the time of matching to classical equations, we may approximate 
the energy density pg- w m^a^. We then promote field cr to a quantum operator: 



&Pcr 



(36) 



The field a can be decomposed into the classical homogeneous background and the quantized perturbation, i.e. 
a = a + 6a. Unlike the inflaton (p whose classical backgroimd is non-zero, because we consider the field & without 
classical background, the leading density perturbation starts with the quadratic rn the operator i5(7^. As with any 
quantum composite operator, we renormalize it with counter terms invariant imder the underlying gauge symmetry 
(here, it is diffeomorphism): 



(37) 



where the subscript r denotes that the operator is a renormalized composite operator, R is the Ricci scalar, and Xi 
are Pauli-Villars fields, which is described in Section |D] We apply this to gauge-invariant isocurvature variable Sg 
defined in Section Hi A I Then we have 



(C) 



(U) 



3H 



3H 



n 



3^. 



(38) 



(39) 



We will not write the hat explicitly from now on. 

In the next subsection, we present how a non-diffeomorphism-invariant estimation of the cross-correlation leads 
to an observationally attractive but grossly incorrect result. In subsections after that, we identify the problems with 
the wrong estunate and calculate the cross-correlation properly. 



A. Plausible but Wrong Estimation of the Cross-Correlation 

In this subsection, we present a plausible estimation of the cross correlation that leads to a large value that is obser- 
vationally interesting. Unfortunately, we will see in later subsections that the estimate presented in this subsection 
can be many orders of magnitude off due to the explicit breaking of diffeomorphism invariance in the treatment of 
the UV physics. Nonetheless, what is presented in this subsection is interesting both as a lesson in field theory and 
as a motivation for the careful correct computation that follows later. 

The isocurvature cross-correlation in the comoving gauge is written as 



« \ ' (40) 
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where we have used dt ( (c^^),,) + 3H ( (c^^),.) ~ for the isocurvature field number density. For an order of mag- 
nitude estimation, we consider a non-derivatively coupled part of the gravitational interaction, 2[,a^SijT'J E Hj^g-a' 
Then the two-point function, shown diagrammatically in Fig. |2j is written in the Fourier space as 



c 



(41) 



xlm 



—3.3^3(^1+^2-?) 
(271)3 



\ki-k2 



l:,{l);i(2) 1 



(42) 



where 



{AB)= / d^xe-'P-'{A{t,x)B{t,0)), 



and Mj; are mode functions for ^ and u, respectively, and 3|'^ means the time derivative with respect to M|:.(f 
This integral is UV divergent, and thus we introduce the horizon scale UV cut-off 



(43) 



10 



(44) 



Moreover, we neglect the contribution from the time range t < tp, where tp is the time when the scale p exits the 
horizon since is oscillatory before the horizon exit. Using the super-horizon approximation for mode functions 
during inflation 



1 H * ,^ . k 



4eMp fc2 



1 f2^-ir(v) 



1 

712 



k 

oH 



1 

712 



2>'+ir(l+v) \aH^ 

where v = -\/9/4 — nP- / H^, the cross-correlation at the end of inflation time is approximately 



(45) 
(46) 



87r2 



TO 

hp 



2 H4 



nit, 



3h2 



(47) 



where we used the relations <C H2 and = H^/AMp£ is the mode function behavior in the long wave- 

length limit. To understand the magnitude of this expression, note that for physical CMB scale comoving momenta, 
we have 



(48) 



for N{p) ~ 0(50). As long as 



1 » ml/H^ > 1/N{p), 



(49) 



It is also helpful to remember that in terms of Fourier space operators/fields, the tilde notation is equivalent to 
where 

A{t,p) = j d^xe-'P-^A(t,x) 

for generic operators/ fields A and B. 
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we can estimate 



-1 

8^ 



(50) 



which is an expression that is valid when the p is far outside of the horizon and a constant H is a good approximation. 
Note that this does not vanish in the limit p 0. We will soon see that this non-vanishing behavior is incorrect and 
is a signal of explicit breaking diffeomorphism rnvariance coming from Eq. ( |44| . Note that if Eq. | |49| is not satisfied 
because nicr = 0, we have 



127t2 



1° 

bp 



In 



-N(p) 



127r2 



(51) 
(52) 



which again does not vanish and is negative. 

As explained around Eq. p4| , the importance of the cross-correlation in the isocurvature bound depends on 
whether /5 is of order 10^^ or larger and not by whether the cross correlation by itself is of the order of curva- 
ture perturbations. To compute /3 defined in Eq. 1 29 1, we need an estimate of {o'^)r correlator which we can take from 
EQ|: 



1 

Ifl^ p 



—f{ma/li,] 



(53) 



where / is a function which can have an exponentially small value owing to the functional behavior 



Combining Eqs. l|29j, l|47|, and (|53|, we find 

^wrong , 



H2 



H 



3 h2 



H 



g3 h2 



(54) 



2 rn^ 
'3 h2 



(55) 
(56) 



which after recalling that N ~ 0(50) and Eq. 1 49 1 gives some hope that a proper computation would give a large 
value for /5 with nicr/H satisfying Eq. |49 1.^^ For example, if |j6| = 0(1), then any appreciable isocurvature perturba- 
tion would be ruled out with the current data, affecting predictions of |'24','30 86 1. 

Recall from Eq. | [24| | that the role of the cross correlation can become important if /3 can become sizable while keep- 
ing a. also sizable. One may worry that the enhancement factor in /5 of Eq. | [53| which is approximately proportional 
to oi may make a negligible in the parameter regime in which j5 is enhanceoTHowever, note that a is controlled not 
just by Eq. ([53) but by 



(57) 



which has a one point function squared in the denominator proportional to the energy density squared of u. One 
can straight forwardly check from Ref. 1301 that the denominator of Eq. | [57| can be tuned such that a can remain 

constant while ((cr^),. (c^^), )p is sufficiently small as to enhance j6 as described in Eq. |56 l 



It is important to keep in mind that we are making an assumption here about the isocurvature evolution when identifying the primordial 
computations of Eqs. j47) and j53| with the CMB observables of Eq. j29| where c, are computed according to the simple transfer treatment of 
AppendixlAl We will discuss this assumption more in detail in subsectionllV A| 
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Given this generic possibility of ruling out a large class of isocurvature perturbation models, we consider below 
the leading gravitational interaction contribution to /5 carefully. We find that unlike the naive estimate given in 
Eq. 1 47 1, there is a suppression in the limit p/{aH) for the mass in the range of Eq. (49 1. The suppression 
in the numerator of j6 precisely cancels the denominator suppression factor coming from jm Eq. (54 1 such that 
no enhancement is obtained, contrary to the naive expectation of Eq. (|56|. This suppression of the numerator in 
the proper computation not seen in the naive estimate can be attributed to a Ward identity associated with the 
diffeomorphism group element of constant scaling of the spatial coordinates. Furthermore, a careful computation 
that we give below will show that the sign of the cross-correlation will be opposite to the naive estimate, owing to 
the fact that the cross correlation here is tied to particle production instead of volume dilution. 

The detailed computation will address also explicitly how same answer to the gauge invariant correlator results 
in two different gauges of comoving gauge and uniform curvature gauge (one can verify this is not obvious from the 
naive estimate presented in this subsection). Another technical care that is taken in the computations below is to ex- 
plicitly specify how diffeomorphism invariant counter terms are introduced to renormalize the composite operators 
intrinsic to Sg. Since the correct answer relies on a gravitational Ward identity, identifying proper diffeomorphism 
invariant regulator and counter terms is important for a trustworthy computation. On the other hand, note that 
the finite parts of the counter terms that remain after the divergences are canceled will not affect the results to the 
leading h expansion that we are concerned with.^^ 



B. Plausible and Correct Estimation Using a Soft-^ Theorem 



Before we describe the actual computation, we give in this subsection a method akin to the soft-^ theorem used 
by Il44l I47H521 |54H60l to estimate the correct answer without a detailed computation. We will also point out what 
ad-hoc ass umpti ons are needed to make this estimate using this theorem. A rigorous computation will be given in 
subsection ( III D I. 

In the soft-^ theorem application to the correlators in inflation, one factorizes N-point function including at least 
one soft external ^ into {N — 1) -point function times the two point function (^^). The well-known example is the 
three-point fimction (^^^) in the squeezed limit in quasi-dS space: 



(2^ 



bp 



2 1 3 

Pa In A: 



TO 

bp 



(58) 



where the superscript on the ^ mode functions denote long wavelength parts. To use this, note that if we neglect 
renormalization of the composite operators, we can write 



d^k2 
(2^ 



72^ 



Using Eq. ([58) and replacing two I, fields with cr fields, we can estimate 



/ 



TO 



d^k2 1 a 
(271)3 fc|ainfc2 



(59) 



(60) 



where the comoving IR cutoff p is required to treat as a constant background field. This effective lower cutoff p 
cannot be justified without explicit computation, but this is physically plausible because {era) does not have any IR 
divergence as long as > 0. One can rewrite the integral in Eq. (60 1 as 



bp 



31nfl 



(f,x) 



(61) 



where the cr^ on the right hand side corresponds to spacetime field (and not its Fourier transform), the p subscript on 
the bracket corresponds to the IR cutoff in the mode function integral, and we assume that there is no contribution 
from the UV cutoff. It is easy to prove that if p ^ is well defined and a UV cutoff is not required, then the right 
hand side of Eq. (|6TJ vanishes in the limit p —> 0. This is in contrast with Eq. (47 1. 



Note that particle production is non-perturbative in h. 
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The vanishing of this function in the p — > limit for > is intuitively understood from the fact that in that 
limit, acts as a spatial diffeomorphism 

f^f(l + ^0) (62) 

(which in turn effectively rescales the scale factor a by a constant factor if we neglect spatial derivatives on long wave- 
lengths) which cannot change (^a^{t,x)) = (^cr^{t,0)) . More explicitly, one can show that the explicit computation 
can be rewritten as 



where 



Q{t) = d^za^{t,)5,jfj{z) (64) 

is the generator of the diffeomorphism associated with Eq. | [62| |. Note that the right hand side formally vanishes 
when the IR cutoff is removed (i.e. p = 0) because in that limit, we find the commutator 

{[Q{t)ro-\t,0)]) ^0. (65) 

This can be interpreted also as a Ward identity. On the flip side, as long as p 7^ 0, ((j^(f, x))^ is not invariant under 

the diffeomorphism Eq. ||62j. The crucial point from this perspective is that diffeomorphism invariance is extremely 
important to see that the cross correlation vanishes for p ^ for a massive scalar field. It is this that one failed to 
preserve in Eq. < |44| . 

As we will show in detail, Eq. | [6T) is consistent with the explicit computation. Note that a couple of assumptions 
that we already mentioned in deriving Eq. |[6T) can only be justified by an explicit computation: namely, the effective 



lower cutoff p in Eq. I6O1 and UV cutoff details associated with renormalizrng the composite operator a^. Such 
complications do not arise in isocurvature scenarios without composite operators. Hence, one of the main technical 
merits of this paper is to provide a explicit justification of Eq. d6l|. Note that because the diffeomorphism gauge 



invariance plays a crucial role in obtaining the correct p dependSice in Eq. (61 1 as explained around Eq. |65}, we 



choose a UV regulator that preserves diffeomorphism invariance in the computation below. 



C. Gauge Invariance of Correlators 

Before we begin our explicit computation, we will check the setup of our computation by demonstrating that 
the manifestly gauge invariant quantities {5gQ and {SgSg) yield the same values in comoving and in the uniform 
curvature gauges. To accomplish this, we use a gravitational Ward identity. 

We first note that the ^ dependent metric perturbations Sg^^"^ and Sg^^"^ differs by a gauge transformation, i.e. 

A^,. = - Sg'fJ = ^^^h) ^ ^ -[£xgh., (66) 



where 



X° = --^, X'=0. (67) 



Their interaction actions differ by 

Their interaction Hamiltonians differ by 

mzaa{t)=H^^l[t)-HfJ^{t) = j d^xa\t)V^'{g,a;t,x)V^X,{t,x) 



(68) 



(69) 
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Then we compare {cr^^y) in the two gauges: 



dt 



f 



dtd^z 



%, V,, {t)T^^ [g, a; t, x)Xy{t, x) 



(70) 



(71) 



where we have integrated by parts and used the quantum version of V^T^'' = 0: i.e. in-ui formalism gravitational 
Ward identities 



iV y{in I Tl'^~^a^ay \ in) 



1 3 

— - z)gf — (m|(7+(7-+|m)^ 



+ 




S (y - z)s,'^}' - — {in\a^ \in) „ 



(72) 
(73) 



whose the notation is explained in Section |B] Note that the remaining term in Eq. | |7T| is a total derivative. Hence, 
we are left with the boundary contribution 



{o-\tf,x)atf,y)f - {(T\tf,x)atf,y))' 



1 



H 



-r2 tOO 



H 



(74) 
(75) 



To make these composite operator correlators well defined while maintaining diffeomorphism invariance (see the 
discussion surrounding Eq. |[65|l), we need a proper covariant regulator, such as the Pauli-Villars (PV) regulator. It is 
straightforward to use the PVregulator here because the above identity holds for PV fields as well. See Appendix [D| 
for a more detailed discussion of the prescription of the PV regulator. 

Using Eq. it is now trivial to show that {SsO^ = i^sO'^ and {3sSs)^ = {Ss^s}'^- Because Ss 3 all {(tI), the 
denominator of this expression also transforms: 



A(5c 3 



Mai 



H la}) 



(76) 



which leads to a cancellation of Eq. 1 75 1 consistently to leading h ^ approximation. Hence, we have a nontrivial 
consistency check of the computationthat we are setting up. 



D. Two-point Functions 



In this subsection, we present a rigorous computation of f> defined in 1 29 1. To this end, we need to calculate the 
two-point function ( {a^)^. C) and ( {a'^)^. (c^^),,) where the renormalized composite operator Il67ll88ti99l is 



[a^)r =[a + Y,Xnf + (A, m^) + 5Z^ (A, m^)R 



{77) 



which is discussed in greater detail in Sec. ID 2 1. Here we are going to use the comovrng gauge for the computa- 
tion because of its advantages that we state below.^'^ As shown in Eqs. 133 1 and 134 1, the gravitational interactions 
in the comoving gauge are derivatively (i.e. p^/a^) suppressed except the (i/) -components. In other words, the 

contributions from T^^Sg^^ and T^'^g^'' interactions are 0{p^/a^), where p is an external 3-momentum. Further- 
more, all counter term contributions are also derivatively suppressed in the comoving gauge: SZq {^) — and 



This computation has been done also in the uniform curvature gauge, which is presented in Appendix |e] Particularly, in the massless limit, 
we explicitly calculate up to the next leading term including all gravitational couplings. This shows that the next leading terms are indeed 
suppressed by the factor p^/fl^. 
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- ^ .C 



P 



ko = p — ki 

Figure 2: Two-point function at one loop order. 



SZi{RQ^ = 0{p /a ). Therefore, we don't need the counter terms to compute the non-derivatively suppressed 
contributions, but we still need a regulator for UV divergences in the computation. The regulator dependences and 
the UV divergences will automatically disappear together in our final result. 

Now we compute the two-point function shown in Fig. |2| which is written in the Fourier space as 



N=0 



O 



(78) 
(79) 



where we have introduced the Pauli-Villars (PV) regulator (see Appendix [D| f or more details) and 



N=0 ^ ^ 



(80) 



where hq and cr„ are the physical field cr and the PV field Xn (here, n e {1, 2, s}), respectively, and s is the number 
of introduced PV fields. 

Interestingly, this integral can be computed in any FRW space-time. We first compute the second term contribution 
in Eq. ( [80) defined as 



Expanding in mode functions, this becomes 



afj{t,x)^{t,0),i^zCN { —ctn 



(81) 



-AC 



N 



d^ki 
(2^ 
ki ■k2 



^Jk25\k^+k2 



dtzClz 



Im 



Cp(0Cp(*z)WN,)ci(0"N,fci(*z)"N,)c2(0"N,;t2(^2 



(82) 



where M]v are the mode functions for fields u^. Because ^ oscillates before and freezes after the horizon exit, we 
neglect the contribution before the horizon exit. Furthermore, we can neglect the 0{p^ / a^) term and factor out of 
the time integral. We thus find 



1(2) (p) 



d% 

(2^ 



^d%S^ih+k2-p) 



X / dt^al ( ^^K^ 



Im 



O 



(83) 
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where tp is the time at which scale p exits the horizon. Note that we drop subscript N and field normalization Cjv for 
convenience, but we will put it back later in the final result. Moreover, we neglect the low momentum phase space, 

i.e. min{ki,k2) < p, because of \u)J^ < 0{k^^) and the spatial gradient factor ki ■ kil cP-. 



ki<v (27r) 



^^d%3Hh +k2- p) f dt.al ( ^ 



Im 



< O 



(84) 



Then the main contribution of the integral comes from the phase space ki,k2 > p, and thus p behaves as an IR cut-off 
(see the importance of this IR cutoff in the discussion surrounding Eq. |[6T|). 
Since ki,k2 > p, we Taylor-expand the integrand with respect to p anatake the leading term. Then we have 



k\ 



Im 



o 



h (2/1)^ h, 

Now we are going to compute the time integral. Recall that the differential equation for mode fimction ui_ is 



(85) 



ill + SHwi- 



Applying to the equation, we obtain 



fc2 



m ] U]^ = 0. 



(86) 



(87) 



where i/j; = g^M/t- Note that the homogeneous solutions for i/j. are uj^ and m|. Thus, we use the Green function 
method to find a solution 



yic(0 = fdt'"^{ui{t)u,{t')-u,{t)uut')) (^-I'pj U,{t'). 



k^' 



From this, we find 



din A: 



|«,(0r = 2ReK(0y,c(0] 



.k\ 



dtzCi^'^lm 





rtp 


f dt^ + 


/ dt. 




J — OO 



4fl:^^Im 



ui{t)uf{t. 



(88) 

(89) 
(90) 

(91) 



The second term is oscillatory with respect to k so that we can safely neglect it after the momentum integral. Inserting 
this back to the integral |[85), we obtain 



-C 



-C 



N 
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2 







2 r d^k d 



{2Tzfdlnk 
" 27r2 



|"N,)c(0|' + o(^^ 



UN,k{i)\ 



3 a: 



'N 



o 



where we have put the subscript N and the field normalization back, and 

f 



'N 



d^k I ,2 

73 \UN,k[t)\ , 



(92) 
(93) 

(94) 



where the subscript p stands for the comoving IR cut-off of momentum. One can then compute the contribution of 
the first term in Eq. dSOl in a similar manner: 
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o 



4(t,f)C(t,0),f 



-3)£,(z)az) 



Hence, we obtain 



N=0 



o 



bp 



27r2 



|"p(0| 



o 



(95) 
(96) 

(97) 
(98) 



where Wp is the mode function for physical field cr. 



Comparing the computation of Eq. | [98| with the estimate in Sec. Ill A we see two crucial differences: 

_-, 2 

1. There is a cancellation of the 3Cj^ 



\ (^n) p ) term that is sensitive to mode summation that extends to sub 

horizon modes. 

2. The Auv dependent term in Eq. | |93) in the prese nt computation disappears after accounting for the PV regu- 



lator fields. In contrast, the estimate in Sec. Ill A leaves behind a A(jy = flH/„y- dependent contribution due to 



the ad hoc nature of the UV cutoff which does not preserve diffeomorphism. 



Finally, putting the results 1 93 1 and d96l together, the two-point function becomes 











TO 











r^(i/)tf 



(2fl(f,)H) 



3-2v 



massive scalar in dS 
massless during quasi-dS 



O 



(99) 



I T-T massies! 

where Hp denote the Hubble scale at which scale p exits the horizon, v = -\/9/4 — nP- /li^, and f ^ reminds us that we 
are evaluating this at the end of inflation. We have applied (quasi)-dS mode function in evaluating 1 99 .^^ One can 
easily check that Eq. ( |47| is consistent with Eq. 160 1. 

As explained near Eq. 1 65 1, the vanishing oflne cross-correlation in the limit p ^ is expected from the dif- 
feomorphism Ward identity. For a nonvanishing p, one might ex pec t the cross-correlation should be 0(p^/fl'^) by 
Taylor-expanding the cross-correlation at p = 0. However, Eq. |99j interestingly shows that the leading term of 
the cross-correlation is not analytic at p = and thus not I fl'^-suppresed. Indeed, for any small p/a{te), we can 
diminish the suppression by making 3 — 2v ^ 0^ through the limit m/ H ^ 0. 

To finish the computation of /5, we also consider the two-point correlator ( (c^^) {cr'^) ^ ) showing up in the denom- 
inator. Again, the comoving gauge is convenient for this computation. Although the correlator is UV divergent, 
because the counter terms associated with the divergence are derivatively suppressed, we do not need to include the 
counter terms in computing the IR contributions and the non-derivative contribution of the correlator is insensitive 
to renormalization. Furthermore, the IR contribution using the super-horizon approximation is not UV divergent. 

That means the UV contribution and the IR contribution are cleanly separated. Thus, we can estimate {{c'^)^ 

using only the super-horizon approximation unlike in the computation of {{o'^)^ C)- We find 



((<72),(C.2)^.)^=2 



d% 
(2^ 



d%s' (/ci + - p) I (0 n (0 r + o ( ^ 



(100) 



After inflation ends at time tg, the cross correlation is expressed as 

((^>P=/r(('^>p 

where fj accounts for the change in the mode-function behavior after the end of inflation. As alluded to in the discussion near Eq. j24| , the 
factor fj cancels out of the expression in fi due to its appearance in the denominator ^A^A^^. The factor fj can also account for the corrections 
in the superhorizon mode function behavior during inflation due to deviations away from the exact dS background. 
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where A/r is a comoving IR cutoff. Evaluating this with dS super horizon modes and assuming m < 3H / 2, we find 
the value at the end of inflation to be 



I All, {2ny 

1 H4 1 / p ^ 

2tz^^3-2v \a{te)H 



2-4+4F |r(i/) 



-21/ 



-2v 



A 



IR 



3-2v 



(101) 
(102) 



In Eq. JTOOl, we have introduced a comoving IR cutoff Ajr which corresponds to the statement that inflationary 
era had a beginning in the finite past. Explicitly, we cannot use the Bunch-Davies vacuum boundary condition for 
modes that left the horizon before the beginning of inflation. This means that 



A 



-(Ntot-Mp)) 



(103) 



where N^ot the total number of efolds of inflation, N(p) is the number of efolds before the end of inflation at which 
the mode p left the horizon: i.e. p/a{N) = H. This cutoff is related to the box cutoff introduced in Il86lll00[|l01l . 
Numerically, A/r ^ p is irrelevant when 



mt, 



H2 Ntot-N(p)- 



(104) 



For situations in which this condition is violated, IR effects are important, and our computation is only qualitatively 
suggestive since A;r has to be resolved using more detailed description of the beginning of inflation. In particular, 
since we do not physically expect N^q^ = oo, nia = situation is not accurately captured by our computation. Of 
course, the IR sensitivity here is not important as far as the importance of the cross correlation is concerned since 
the qualitative behavior of ha ving p/A/R oo is to make the correlation even larger making the fi parameter even 
smaller. Finally, note that Eq. 1 104| can easily be more stringent than Eq. | |49| |. 
Hence, we conclude 



nig 
" H 



In-J^ 



-1/2 



massive scalar in dS 



massless during quasi-dS 



(105) 



where for the massive scalar case is assume to satisfy Eq. p9) . Although this in principle is a generic prediction of 
isocurvature scenario, the magnitude of around 10^^ is difficult to probe experimentally since the current sensitivity 
is at the level of lO"^. 



IV. APPLICATION 



The /3 computation presented in Eq. 1 105 1 is not sensitive to pa that is involved in the definition of the isocurvature 
perturbation 5^ . Instead, it is a property of quadratic nature of the scalar composite operator during inflation. Since 
Eq. | |105| does depend on the masses, in this section, we motivate couple of the mass parameters from well-motivated 
nonthermal dark matter models: WIMPZILLAs 1102111081 and axions f 109-1 11]. Although these two particles have 
different physical origins, they share some common properties as a cosmological component. Firstly, since they are 
massive (at the CMB time at least) and weakly interacting, they both are good CDM candidates. Also, they can be 
gravitationally produced during or after inflation, and this gives rise to isocuvature from their density perturbations. 
Furthermore, when their background field values are negligibly small, the isocurvature perturbation from these 
particles is approximated by quadratic form a^. In that case, they would present detectable non-Gaussianties | 
|S6 1 and their cross correlation is characterized by Eq. ||105|. 



A. Weakness of a Interactions with ip 



To connect our computation of j6 to observables, a post inflationary isocurvature scenario is necessary. For the 
illustrative situations of axions and WIMPZILLAs, it is sufficient to assume that a has an extremely weak interaction 
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with the reheating degrees of freedom xp and the rnflaton (p such that the transfer function of a is trivial after inflation: 
with sufficiently small interactions, a and /3 of Eqs. | [28) and ( p9| computed during inflation can be directly matched 
without any further transfer function computations to isocurvature initial condition for CMB codes such as CMB- 
FAST. In this section, we quantify the requisite weakness of the interactions and qualitatively discuss the situation 
when the weakness assumption is invalid. For example, we will show below that ordinary WIMPs are too strongly 
interacting with the reheating degrees of freedom for this assumption to be valid while axions and WIMPZILLAs 
are sufficiently weakly interacting. We also qualitatively describe what extra work needs to be done to apply this 
paper for observations in situations in which the dark matter particles are not extremely weakly interacting.^^ 

At the linearized classical equation of motion level, we have the gauge invariant perturbations {Q} being gov- 
erned by a linear time evolution operator 

0[{Q}] = (106) 
where the initial condition for the isocurvature species / = is given by 

Utd=fiti) (107) 

which in turn is set by the inflationary physics. For example, the initial time f, can be set to be the time of end of 
inflation. The final Co-{tf) will contain contribution which does not vanish in the limit / — )• 0. Hence, one can write 

Utf) = Gf^lfitdA + Gr^[0,C;^,(t,)] (108) 

where Gf [D] is the cr component of the Green's function derived from the linear operator O which takes the initial 
data D and maps it to the final value of ^cr{tf)- Note that we have implicitly assumed the boundary condition such 
that Gf^ [0, 0] = which means that Gf^ [/(i,), 0] vanishes as / (f,) 0. 

Now, we will consider two situations in which boimd the picture of super weakly interacting scenarios. In the first 
scenario, the thermal plasma generated by the inflaton decay will interact with cr sufficiently strongly to make mix 
strongly with ^. In the second scenario, the inflaton decay to a directly will realign a fluctuations during radiation 
domination to those of ^, even though u and reheating products are not interacting appreciably. 

First, consider the effects of radiation dominated thermal plasma on a. The mixing rate governing Gj^ [0, ^y^cr(f;)] is 

the production rate of a particles from the thermal plasma. Typically a single channel involving particle y dominates 
the production of the a particle from the plasma. (If there are more channels, the discussion below can easily be 
generalized.) We thus expect a qualitative behavior of 

Gf^[0,q^,{U)] ~ (^1+tanh 

where T{yy acr, trnax) is the reaction rate for this process at the time that the production rate is maximum (in 
T{yy — > aa, t) is maximum ait = tmax where fmax G [f// f/])/ is the expansion rate, and = O(^tot)- 
Hence, one sees that the information about the isocurvature perturbations depend not only on 

T{yy cra,tmax)/H{trmLx) 

but on f ^ since frnax is restricted to be in the range frnax G [U>tf]- For example, the usual CMB code is run starting 
with an initial condition at T ^ ^BBN' This means that f/ ^ '^BBN required to use the inflationary correlator 
computations in the CMB code. A general computation of needed for the prediction of isocurvature perturbation 
effect on CMB temperature is beyond the scope of current work. To be able to trust the trivial transfer function of 

Gl[f{ti),0] « f{ti) » Gf [0,q^^{ti)], (110) 



r(yy aa,tmeLx) 
H{traax) 



(109) 



Because of the cross correlation result in this paper is small, the discussion here is a bit academic if this discussion applied only to the cross 
correlation result. However, the discussion here applies to the isocurvature 2-point function found in the literature (9 30 41-43 100llll2l|113l 
which has a realistic chance of being observable in near future experiments. 

In our scenario, the isocurvature species stand for the degrees of freedom constrast with the radiation degrees of freedom. 
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for superhorizon modes (where ti is say at the end of inflation^''), we can require 

Tjyy ^ cr(r,tmax) ^ CaiU) ^^^^^ 
H(fmax) ^tot 

where tmax can be at any time between inflation and the time at which botmdary conditions are set for the CMB 
code. This sets a bound on the cross section {av) for yy aa to be 

{(TV) « /"^['V*^' ' 4.2 X 10-25 GeV^2 (112) 



^ioiih) gy VlO^GeV 

where the bound becomes more stringent for higher reheating temperatures. 

This number should be compared to typical thermal WIMP DM candidate annihilation cross section of 10^^ GeV"^ 
and a high energy s-channel scattering at Trh mediated through a vector boson with a dimensionless coupling 

Oil 

(t^^yy^A/' flight states) ~ ^ (113) 



Hence, one sees that WIMP dark matter cannot play the role of the isocurvature perturbations. That is why if we are 
to identify our computation of a and /5 directly to physical observables, we have to choose the isocurvature degree 
of freedom to be nonthermal. 

Even though the current work applies most immediately without changes to nonthermal dark matter scenarios 
having extremely weak interactions, Eq. | |112| is still much bigger than gravity mediated s-channel interactions 

(-y.-..-^^> - i^lf (115) 

For example, axion cross sections for gluon coannihilation behave as II115I1 

(-«.^x> - (117) 



1012 GeV 



2 

10-^» GeV-^ (118) 



/" \ in— 28 r^„\/— 2 



where /« is the PQ breaking VEV. Hence, there is a large class of weakly interacting models for which this work di- 
rectly applies without modification. For models for which Eq. | |112| is not satisfied, one needs to compute the transfer 
function associated with the mixing. Nonetheless, this work will still be useful in setting the initial conditions for 
such computations. 

Let's see qualitatively what happens when Eq. 1 112 1 is not satisfied. In that case, we expect mixing between 
isocurvature and curvature perturbations 

Uif) = Gr,[/(i,).0] + G^;^[0,^j^^{ti)] ~ 0(^.) +0(^tot)- (119) 

Since the curvature perturbations will analogously be 

aif) = Gf^\f{h)A + Gf^[G,Zj^,{U)], (120) 



Note that as discussed in footnotelwl a can also receive corrections from the departures from the ideal dS mode function evolution as well as 
from the time when m/H becomes larger than unity. As discussed there, the quantity fi is not as sensitive to these corrections. 
Similar arguments can also be made from unitarity 11141 . 
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we would then have 



= 3 {{Gf^[f{U),0] - Gf^[f{U),0]} + {g^^[0,Cj^,{U)] - G« [0,^^-^,(f,)]}) 



(121) 



Up to the accuracy that all species are equipartitioned, this quantity may vanish since there is cancellation in each of 
the terms in the bracket. It is beyond the sc ope of the current work to compute more precisely this cancellation we 
are focusing on scenarios which satisfy Eq. 1 112 1. 

Suppose there is a direct decay of the inflaton to a, and suppose there is no other appreciable interaction between 
u and other decay products of the inflaton. In that case, it is better to set the initial time f , to be at the time of inflaton 
decay completion such that Gf [f{tj),0] is still trivial. In that case, we have 



+ 



^pfe'-^v) ^ ^^(decay) 



(grav) (grav) 



where r/ has been defined in Eq. (TT) . Hence, we have 



(decay) (decay) 

Ccr 



(decay) -(decay) _ . 



If ^^^^^^^ = is assumed, then 



(decay) 



(grav) 



(122) 
(123) 

(124) 
(125) 

(126) 



This equation says that if most of the inflaton energy density goes to cr, then the isocurvature is negligible. 

In the next two subsections, we now consider couple of mass motivations for nonthermal dark matter isocurvature 
candidates. 



B. WIMPZILLA 



The WIMPZILLA was originally proposed to avoid the restriction from the assumption that the dark matter is 
a thermal relic. Thus, the WIMPZILLA is supposed to either be very heavy and/or very weakly interacting. In 
particular, we consider the possibility that the WIMPZILLA is gravitationally produced during the phase transition 
out of the quasi-de-Sitter phase of inflation. In that case, the model is controlled by two parameters: the ratio of mass 
to the Hubble scale of inflation otx / Hj^f, and the reheating temperature Tr^, where X denotes a massive scalar field. 
Since the energy density is approximated as ~ m^X^ the relic density of X is estimated as 



where we have assumed that ntx ^ He, because a priori we know that we can find proper isocurvature and relic 
density in this mass range. (For a more detailed discussion of the relic abundance, see for example [30].) The 
isocurvature power spectrum depends on the details of the evolution of the background during inflation because 



a . , in-l 



the mode function of massive particle decays as a"^^'^^' (see a related discussion in footnote 17 1. However, we can 
generally obtain a ^ 0.067 if mx ^ Hinf, where H,>,| is the Hubble expansion rate when me CMB scale crosses 
the horizon |30|. The WIMPZILLA isocurvature has also the quadratic form like the axion. It thus generates the 



observable non-Gaussianities estimated as Eq. dSTb. Eq. 1 105 1 translates to the fractional cross-correlation of 



3wjmpz;lla « -0.4^, \/A^ (128) 



which justifies the constraint used in fSO]. Since the naive estimate of Eq. | [56) gives a gross overestimate j6, one of 
the merits of this paper is to put such worries to rest through the proper computation. 
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C. Axion 



In this subsection we assess the relevance of Eq. | |105[ l to the axion scenario. Firstly, we review the axion scenario. 
In 1997, Peccei and Quinn proposed the global U(1)pq symmetry in order to solve to the strong CP problem in 
the QCD. The axion is the Nambu-Goldstone boson associated with the symmetry after it is broken spontaneously. 
Many mechanisms have been proposed to produce axions in the early universe. We focus only on the "vacuum 
misalignment" mechanism here following Refs. |43j Ill3i llT6l4l2 2l. In early universe, the axions are effectively 
massless and gain their mass when the QCD anomaly term (which explicitly breaks PQ symmetry) becomes physical 
after the chiral symmetry breaking QCD phase transition. After the universe cools down and the Hubble friction 
drops below the axion mass, the axions begin to coherently oscillate and they contributes to the CDM component of 
the universe because of their long lifetime. 

Let us denote the PQ symmetry breaking scale by fa- Because o< 0-^ where 6 is the axion angle, the relic axion 
density is estimated as 



2xlO'(ioi^)'^'(^') for Tosc >Aqcd 
5 X 10^ 



fa/N \ 



(02) for Tosc <Aqcd, 



(129) 



where we have neglected 0(1) factors due to diffusion, anharmonic correction, and temperature-dependent mass 
correction, and Tosc is the temperature at which the axion starts to oscillate. The axion isocurvature in comoving 
gauge is written as 

92 -(02) 20i39 + Se^ 



AC) 



COn 



COn 



(130) 



where coa = n„ / CicDM, is the average of initial QCD vacuum angle 6 over the observable universe, and S6 is 
inhomogeneity of 6, i.e. 6{t, x) = 0,(f) + S6{t,x). Then the isocurvature power spectrum becomes 



3.5 X 10 



10 



' fn/N \ 

lois GeV J 



7/3 



where 



A01 



F for/a/N>6x lO^^GeV 
for/a/N<6x lO^^GeV, 



(131) 



(132) 



Since our primary interest is in the cross correlation with 6, ~ 0, we set it to zero. 
Therefore, the adiabaticity parameter a defined in Eq. |[28) is estimated as 



a ^ cv, 



,I1,3x10»(^)"'a 



8 1 X 1017 ( fa/N Y .2 



for fa/N > 6 xlO^^GeV 
for/a/N <6 X 10^^ GeV, 



(133) 



27r2 

k 

N 




(134) 



where Hp is the Hubble scale at the horizon exit of mode p, and A/r is an IR cut-off. Here we have used Eq. 1 102 1 with 
the assumption that the axion is effectively massless during inflation. In the case that 0, <i; dO, the isocurvature has 
the quadratic form of gaussian variable SO, and it naturally becomes non-Gaussian perturbation. The isocurvature 
non-Gaussianity is estimated as Eq. | [3T| . 

These parameter constraints and predictions 1 129[ , 1 133[ and | [3T| already have been investigated in the literature 
I^T, 43 "861 ITT3I with the assumption that the axion isocurvature and the curvature is uncorrelated. Our result from 
Eq. |[105) is 



In 



A 



-111 



< 2.5 X 10"^ 



(135) 



which is consistent with the assumptions made in the literature. 
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V. SUMMARY 



In this paper, we have presented the first explicit computation of the gravitational interaction contribution to 
the cross-correlation between the curvature and quadratic isocurvature perturbations (which include dark matter 
isocurvature candidates such as axions and WIMPZILLAs). Since the necessary and sufficient condition for the 
cross-correlation to dominate over the isocurvature perturbations in the temperature two-point function is > 
4 X 10 we have explicitly computed /5, which incidentally is not sensitive to the background number density of 
the isocurvature degrees of freedom and post-inflationary mode function changes on superhorizon scales. Although 
a naive estimate of /5 based on a diffeomorphism violating UV cutoff leads to the possibility of /5 0(1) due to 
a large ratio that can appear between the numerator and the denominator of the expression for j6, our explicitly 
diffeomorphism invariant computation leads to |/5| < A^/2 ~ 2.5 x 10^^ because the numerator has a suppression 
as a consequence of a diffeomorphism Ward identity. Unfortunately, this is far below the current observational 
sensitivity of |,6| > lO'^. 

The smallness of the cross-correlation is explained by the fact that the super-horizon mode of the curvature per- 
turbation ^ can be smoothly connected to the gauge mode, which is the spatial dilatation, in the zero external mo- 
mentum limit. Hence, Eq. 198} vanishes when p = and m 7^ 0. In other words, this can be seen as a suppression 
due to a diffeomorphism Ward identity (i.e. uniform spatial rescaling invariance). A nontrivial structure revealed 
through our explicit computation is the suppression's non-analytic structure with respect to p: the cross correlation 
cannot be Taylor-expanded at p = 0, and this contribution is not p^/ a^-suppressed. 

Our rigorous result which incorporates UV renormalization of the composite operator in the curved background 
is also shown to be consistent with an estimate based on a soft-^ theorem, which allows one to factorize (^^) from 
(y^C) explained in Eq. 161 1. However, Eq. |[6T| re quires two assumptions that can only be justified by an honest 



computation such as what is presented in subsection III D 



1. There is an effective IR cutoff of p in evaluating (u^) due to the external momentum p inserted into the com- 
posite operator. 

2. The only UV renormalization property of ((J^) that is relevant to leading approximation is the preservation 
of diffeomorphism invariance. 

Note that the proper diffeomorphism invariant UV treatment also allowed us to demonstrate that the cross- 
correlation is indeed gauge-invariant with one-loop correction through the gravitational coupling. This gauge invari- 
ance is checked explicitly by computing our cross correlation in both the comoving gauge and the uniform curvature 
gauge. 

physically, the curvature perturbation ^ can affect the particle density per and generate correlations only at its 
horizon crossing, because ^ freezes out after its horizon exit, after which it can be effectively treated as a gauge mode. 
Positive cross correlation corresponds to the situation in which the 1 + ^ enhancement in the expansion enhances the 
particle production (assuming that this enhances inhomogeneity) while the negative cross correlation corresponds to 
the situation in which the 1 + ^ enhancement in the expansion dilutes the particle inhomogeneity. The latter dilution 
effect leads to j6 > 0, while the particle production enhancement effect corresponds to the quadratic scenario that we 
were interested in this paper. This explains the sign /5 < of our result. 

Given the robustness of the smallness of /5, the gravitational interaction contribution to the cross correlation should 
be negligible in most nonthermal dark matter isocurvature scenarios. In addition to giving a concrete computation 
that supports this statement, our work serves as an interesting lesson in computing correlators of composite opera- 
tors in curved spacetime in the context of inflationary cosmology. 
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Appendix A: Behaviors of Transfer functions for Adiabatic and Isocurvature initial condition 

The CMB temperature fluctuation with the leading order approximation (the integrated Sachs-Wolfe term is ne- 
glected) in the Newtonian gauge {B = F = 0, E = 2^, A = -2Y) is 

AT 1 

-Y^lM+^l' (Al) 

where the perturbations on the rhs are evaluated at the recombination. We can obtain these perturbations by solving 
the Einstein and Boltzmann equations with given initial conditions. A projection from a given initial condition to 
the final CMB temperature fluctuation is called transfer function. In the following subsections, we calculate that the 
^-dependence of the transfer functions for the adiabatic and the isocurvature initial conditions. In particular, we 
show that the isocurvature transfer function has the additional suppression factor kgcj I k compared to the adiabatic 
one for small scale k/.^. Here we basically follow the calculation by Ref . II123H124I . 



1. Perturbation Equations 

For explicit computation, we choose the Newtonian gauge for the scalar metric perturbation For simplicity, we 
consider only photon and CDM fluids, which are denoted in the following equations by subscript 7 and m, respec- 
tively. This assumption is valid for the sake of identifying the difference between transfer functions for adiabatic and 
isocurvature initial conditions, although baryon and neutrino should be taken into account for accurate description 
for transfer functions. 

The conservation equations for dark matter and photon fluids in Fourier space are 



6' 



k^Vm + 3Y, 



4 

1 



-A:2yT,-h4Y', 
^7 - 



(A2) 
(A3) 

(A4) 
(A5) 



where ' denotes the time derivative with respect to conformal time r],^ = a' / a, 5a = 5pal Pa- Note that O = Y since 
they are perfect fluids. Vx is the peculiar velocity for fluid X. These four equation are combined by eliminating Vx/ 
and we have 



(fl(^;,-3cD')) = ak^^, 



^11 



40" - ^ (S, 



-4$) . 



The evolution of the metric perturbation is encoded in the Einstein equations. (00) and (ii) components are 



2Mj 



{pmSm + Pj^y) , 



1 2 . 

a pjOj, 



6Mj 



Combining with other components, we also find the Poisson equation 



(A6) 
(A7) 

(A8) 
(A9) 



^7 "7 



With the definition of isocurvature | [20) in Section [Til 

Ss = S„ 



(AlO) 



(All) 
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where we have used pj = Pj/3 and pm = 0, we rewrite the differential equations of fluid and metric perturbations 
in terms of O and 



where 



+ 3'h(i+ c^) CD' + [27^' + 3c2) 



y 



V2 



CD + fc^c^cD 



2 c2 2 . 



3 ( 1 + 



(A12) 
(A13) 

(A14) 



In — )• limit, Eqs. I A12 1 and I A13 1 admit two linearly independent solutions cI)(A:, // — >• 0) = cD' {k), {k, ^ — )• 0) = 
0, and ^{k,ri — > 0) = 0, 3s{k,r] — )• 0) — Sg{k), which corresponds to adiabatic initial condition and isocurvature 
initial condition, respectively. 



2. Adiabatic Initial Condition 



For large scale perturbations, which enters the horizon later than the recombination. remains zero according to 
Eq. | |AT3) , and thus Eq. | [AT2) is rewritten as 



rf^cD 21y2+54y + 32 d<i> 



CD 



0, 



dy^ 2y(y + l)(3y+4) rfy y(y + l)(3y + 4) 
where is called as Kodama-Sasaki equation. This differential equation can be exactly solved, and we find 



(A15) 



^ikuy:$>l) = j^^Hki), 



(A16) 



where the subscript I stands for "super-horizon". For photon energy density S^, Eq. |A4l in the long wavelength 
limit yields 



and also Eq. dASj gives 



1 

-Sj — cD = const. 



Syiki^ij ^ 0) = -2^iki,tj) = -2^'iki). 



(A17) 



(A18) 



For small scale pertu rbation, which enter the horizon during the radiation dominated(RD) era, in the early RD 
limit t] "C t]eq, Eq. |A12l becomes 



4 k^ 
CD" + -CD' + —CD = 0, 
V 3 



(A19) 



and its solution with the adiabatic initial condition 



{sin wrj — wf] cos wtj) cD'(fcs), 



where w = k/ -\/3. After the perturbation enters the horizon, 

>f.n \ 3cosii;j; , 



[wf]) 



2M^ 

l^{ks,ri) = 6^' (ks) cos ivT], 



(A20) 

(A21) 
(A22) 
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where the subscript s means "sub-horizon", and the second equation is obtained by the Poisson equation < A10| . 
Plugging this solution into Eq. ||A6|, we find that 



(A23) 



Sm{ks,^ <Veq) ^ -9*'(fcs) [inWf] + J 



2 ' 



where 7 is the Euler Gamma constant. This shows that the dark matter density perturbation grows logarithmically 
during the RD era. 

Now we should match this with the solutions in the matter dominated(MD) era. Because the time derivatives of 
O is negligible compared to the spatial derivatives, Eq. I A6 1 is approximated as 



where we have used the Poisson equation I A10| . Then, it is rewritten as 



y(i + y) 



and its general solution is 



i + 2y 



3 \ dSn 



1 + oy In 



dy 2 

y/i + y + i 
v/r+y-1 



0, 



(A24) 



(A25) 



(A26) 



Matching this solution with Eq. I A23 1 at y -C 1, we find 
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y^'{k) ( ln2w;;;* + 7 - - ) wheny > 1. 



(A27) 



where = rjeq/ (^v 2 — 1 j = 2Teq. Note that we have used the results from the Friedman equation 



y 



y y^ 



(A28) 
(A29) 



and Eq. \A27\ corr espon ds to Eq. ( 150) in Ref. IIT24ll . 
Then using Eqs. I A10| and I A27| , we get 



■».(/c.,,>,„)=. '"'°f-";'^ -fe). 

{0.27ksrjec,) 



(A30) 



This shows that the gravitational potential is frozen after the matter-radiation equality. Similarly, we first find the 
general solution of Eq. jATl for sub-horizon modes 



ci cos wrj + C2 sin wrj — 40, 



(A31) 



where we have neglected that time derivatives of <$>. Then matching this with Eq. | |A22| , we get 



6^{ks,r] > r]ea) ^ 6 cos [ivt]) - A— ^ -f 

{G.27ksr]ec,) 

Now we return factors due to the Silk damping and the acoustic sound speed 



^\ks). 



(A32) 
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4 In (O.lSfcs?.^) 
3c| (0.27fc,;?,,)^ 



1>'(/c«), 



(A33) 



which is Eq. (153) in Ref. II124I . Notice that the the first term is dominant for the scales we are interested in. However, 
the second term becomes important for very small scales where the diffusion damping is not negligible, k^ku- 
Finally, the SW term dAll becomes 



AT 



6<i>' {k) cos zvrj iik > kcq 
Acl>'(fc) if 



(A34) 



Note that 



(A35) 



3. Isocurvature initial condition 

For large scale perturbations, 6$ remains constant, and Eq. I A12| has the solution 

cD(A:,,//) = 



x\ + 6x + 10 „• , 



{x + iY 



where x = t] / rii,q.lr\ the MD era, Eq. I A36 1 gives 

1 



&m{ki,r] » rjeq) = jSj{ki,r] » j]^^,) 



mi), 



(A36) 



(A37) 



where the last two equations are obtained from Eq. | |A8| . 

Now, we will see how the perturbations evolve during the RD era, and how they are connected small scale pertur- 
bations. In the early RD era, the source term and the last term on the left hand side of Eq. I A13[ is negligible because 
they are high er ord er in y. Thus, the solution 6$ remains constant even inside the horizon. In tha t case , Eq. |A12l 
becomes Eq. I A19 1 with the source term Sg/lyrj-^^. Then we find its solution that matches with Eq. I A36 



V 



Furthermore, in the wtj ^ limit, we have 



1 + 



{zvt]) 



(cos WT] + wr] sin a;// ) 



(A38) 



and putting this into Eq. ||A8|, we find that 



18 



7 



(A39) 

(A40) 
(A41) 



As explained in Section II B we have that O and 5^ grows like a during the RD era, meanwhile 6m decreases. 
For sub-horizon modes, Eq. I A38 1 becomes 



y 

3 V 2 



sirvwrj j 6^g{ks), 



(A42) 
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and again plugging this into Eq. I A8 1 yields 



3 sin wt] 
2 wrj 



wt] 



(A43) 
(A44) 



Matching these with general solutions of perturbations I A26| and I A31| , and also using Poisson equation I A10[ in 
the MD era, we get 



Then the SW term becomes 

AT 



1 



0.35ksf]eq 
1 

{0.8k,t]eq) 



■ sin {wt]) +4 =■ 

{0.8kst]eq) 



-^^^5i.{k)sm{wt]) ifk>keq 

-l3'sik) if k<t]-\ 



(A45) 
(A46) 
(A47) 

(A48) 



Now we see from Eqs. I A34l and I A48| that the isocurvature transfer fimction has the additional suppression factor 
fcgg / A: compared to the adiabatic one for small scale k > k^q. 



Appendix B: Review of Diffeomorphism Invariance 

A symmetry in a classical field theory is preserved at the quantum level, if the regulator preserves this symmetry 
and if the functional measure is invariant under the symmetry transformation. The quantum symmetry is reflected 
in the transformation of the correlation functions. 

For example, consider a scalar field cr on a fixed manifold {A4,g). The two point function is 

{<7ix)a{y))g = f D<pe''('^-'S)cr{x)a{y) (Bl) 



The two point function only depends on the metric field g and points x, y. Intuitively, the symmetry says for any 
diffeomorphism cp : A4 ^ M., the metric field and the points changes as 

g g = (<p-i)*g, x^x = (p{x), y^y = (p{y) (B2) 
then the two-point function should remain invariant, i.e. 

{cr{x)a{y))g = {a{x)a{y))g. (B3) 

The Ward identity is the infinitesimal version of this relation. 
Let (p = exp(eX), then 

g = exp(-eXrg = g-e£xg+--- (B4) 



5{g,cr) = S{g,a)-e j d'^x^'^TrCx{g)iiv + ' 



(B5) 



cr(x) = (7(x) + e£x(7(x) H (B6) 

Plugging this into Eq. | |B3[ and Taylor expand with respect to e, one get 

- i j d^Zy^^£.x{g)iiv{z){Tt'crxCry)g + {Cx{o-)xO-y)g + {crjCx{a)y)g = 0. (B7) 
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Or equivalently, using 

>Cx(^)^v = V;,X, + V^X^ (B8) 

and perform integration by part, we obtain 



which is the Ward identity for the path ordered vacuum expectation value. We can then write down the in-in 
expectation value Ward identity as 

iV^{in\Tr^a+a+\in)g = -^5^{x - z)g'^'' ^{in\o-+(T+\in)g 

+ -^^^(y-z)g-^(m|<7+^+|m)^ (BIO) 

iV^{in\Tr~cr+a+\in)g = (Bll) 

where we kept the external operator inserted on the forward branch. The fact that Eq. | |B11| has no contact term is 
easy to understand, since ^ is inserted on the backward time branch of the manifold, it can never contact points 
X and y. 

Appendix C: ADM formalism and Interaction Hamiltonian 

We consider an inflationary model with the inflaton (p and an extra free massive scalar a, where cr is only gravita- 
tionally coupled with cp. 

S = I {dx)^-MlR+[-^g^'''^^,ct>^y(l> - Vicp)] + [-^-g''''d,,CTd,a-Uia)] (CI) 



where Mp = = 1 and (dx) = d^xy | det(g^„/) | . The metric can be parametrized using ADM formalism 111251 ^^, 

^f'^ { h.jNi hij )' ^ \ N'N-^ h'i - N'N'N^^ J ' ^ ' 

where hjj is the metric tensor on the constant time hyper-surface, and h^l is the inverse metric. We use Latin indices 
i,j ■ ■ ■ for objects on the 3-dimensional constant time hyper-surface, and we use hjj and h'l to raise and lower the 
indices. Then the action | |C1) is rewritten as 

S = ^J{dx)Vh NR^^^ -2NV{(p)-2NU{cr)+N-^ (^EijE'' -E^^ +N^^ (j>-N'diipy -Nh'idiCpdjCp (C3) 
(cr - N'dicr^ ^ - Nh'idicrdjii , 



where E/y and £ are given by 

-7(3)m._V7(3) 



Erj = ^{h.j-Vr^j-^j^')- (C4) 
E = Eijh'i. (C5) 



Consider the background solution driven by the inflaton. 



=*('). = ( 4,,. ) . (C6) 



We use ( — h ++) sign convention for the metric, and physical time t . 
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where they satisfy the background equations of motion 

H = - 
^ + 3H$ + V'{4>) = 0. 



(C7) 

(C8) 
(C9) 



The action for the perturbations can be obtained by Taylor-expanding the full action around the background solution. 
However, we may reduce the number of variables by imposing the ADM constraints: 



= l[R(3)-i^(E^.^.E'7_E2)]_2NT00 
N ' N 



where 



JflV 


rj,}lV rj,}lV 








= -^''^ 








la 






hU(cr) 





(CIO) 
(Cll) 

(C12) 
(C13) 

(C14) 



and choose a gauge. 

One commonly used gauge is the comoving gauge, defined by 



where 



5cp = Q, 7;,. = 0, 9/7,7 = 



The solution of N and N' is 

We find the scalar metric perturbations are 



-2i 



4 +^#10. 



(C15) 



(C16) 



(C17) 



(C18) 



where e = H/H^. Plug ging in the linear metric perturbation back to the action JC3l, we can get the perturbed action 
action up to cubic order 



S(C) 



AC) , AC) , AC) , AC) , AC) , 
J-yj T '-'^^^ ~r ^aa ' 



where 



AC) 
AC) 

^laa 



(C19) 

(C20) 
(C21) 



In this section, Latin indices i,j are raised and lowered by i/y, and repeated indices are contracted. 
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The ^ cubic interaction and graviton actions can be found in 1*441 . 

Another commonly used gauge is the uniform curvature gauge, in which 

h,j = a^t)[e'^]^., ju=0, 3/7/; =0. (C22) 

In this gauge, the inflaton degree of freedom is in Sep. However, this degree of freedom can be represented using the 
gauge-invariant variable 

^=-^3(pW (C23) 

In this gauge, the ADM constraint renders 

N(i'^) = -eC, Nf'^'=a,[e^e] (C24) 

We get the linear metric perturbation as 



The free action is the same as in Eq.| C20| , and u-^ cubic interaction action is 



4S = fd'xaliT^'e^+T'^'e^a (C26) 

From these perturbed actions, we can obtain the interaction Hamiltonian. Particularly, note that up to the cubic 
interaction, £,„f = —Tiint- Thus Si^^g- = — J dt H[^^^{t). 

Appendix D: Renormalization of Composite Operators 

In renormalized perturbation theory, one requires a regulator and renormalization condition. In order to preserve 
the diffeomorphism invariance, we need to adopt a covariant regulator. Here we choose Pauli-Villars (PV) regulator. 



following 1 126 . 127J. We will first review PV regularization in subsection D 1 and renormalize in subsection 



For correlators involving time integrals, we describes the adiabatic expansion of time integral in subsection I D 3 1. 



D2 



1. Pauli-Villars Regularization 

We introduce a set of scalar regulator fields Xn for n = 1, ■ ■ ■ ,s with the following free Lagrangian 

r-PV = E C„ (-Ig^-'dyXndvXn - \Mlxi) ■ (Dl) 

11=1 \ ' 

The number of regulator fields s depends on how many independent divergences one need to remove. In order to 
eliminate UV divergences up to some even order 2D,we must take the C„ and regulator masses Mn to satisfy 

EC^'^O' E C^,1m2, = 0,--- (D2) 

N=0 N=0 

EC-^Mf = -mf (D3) 

n 

where we used the notation Mq = and Co = 1, and let = a and tJ,, = Xn- We use A to represent the set of Mn, 
and the regulator dependence should be removed by counter terms when M,, goes to oo together. 
On a homogeneous FRW background, the physical and regulator scalar field can be quantized as 

WN,irM\ = ia~^{t)3^{x-y)SNMC^^ (D4) 
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with the following mode decomposition 



crN{x,t) 



where U]^^p{t) satisfies the usual equation of motion 



(271) 



fc2 



UN,k + 5Huf^^k +\-^+^n] ^N,k = 



with the Bunch-Davies initial condition 
and Wronskian conditions'^ 



— exp ( —i f —^dt' ] for f ^ — oo 
2ka{t) ^\ J a{t') 



^N,k^*N,k ~ ^T^,k^*N,k — • 

Because M„ ^ H, Eq.| |D7| possesses the WKB-type solution 

un,k{t) = . ]. 3 , . exp (-i fcv,{t')dA 1 + + Att^ 



(D5) 
(D6) 

(D7) 

(D8) 

(D9) 

(DIO) 



where coj^ = \/k?- 1 cP- + and /, are of zeroth order in a;j-. Since we have to regulate up to quadratic divergence in 
correlator computations, we need to know 



|wia(0| 



^^2Re^^|/,(0|V2ReM)^ 3) 



<^k{i) 



up to second order. Due to the equation of motion | |D7| , fi should satisfy 

7i 



(Dll) 



Also, the Wronskian condition ||D9| yields 



L (ii I 2H^ I + 5H'M^ 

2 a;? 4: cof 



Then plugging these two results to Eq.| |Dll^ gives 



\u„,k\ 



2cO]^a^ 



H + 2H2 (H + 3H2) m2 SH^Mf! 



■0(lo- 



(D12) 



(D13) 



(D14) 



2. Renormalization of Composite Operator 

The renormalization of composite operators in curved space-time is the same as in flat space-time(see e.g. f88','89j 
11271 ) , just with new possible counter-terms made from curvature tensor. For an operator of dimension n, one need 



Our treatment here differs from [1261 in that the physical scalar field cp here has no background solution, and the regulator field Xn does not 
mix with (p by mass term. 
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to consider all possible counter-terms of dimension n or less. In our example model with free massive scalar a, we 
renormalize cr^ as 



= {a + J2Xn)^+SZo{A,m^)+SZ^{A,mcr)R, 



(D15) 



where R is the Ricci scalar 

Next, we compute SZ/s divergent part. For example, let us consider the one point function 



(D16) 



In order to determine the counter terms SZq and SZ^, we introduce a comoving scale Q such that H ^ Q/a ^ M„ to 
break the Fourier space into the UV and the IR sector Then we use the WKB solution | D14| for A: 3> Q. Furthermore, 
the contribution from the PV fields for A: ^ Q vanishes since it is suppressed by 1 / Mn ■ 



N 



Q d\ , 



N=0 



(271)3 



^"i' dh . 1 2 
a 10 



9671- 



N=0 



N=0 



Note that the arbitrary comoving scale Q in the first two terms should cancel each other. 
In order to absorb the PV regulator dependence, we need 



SZo 



SZi 



1 



167r2 
1 

967? 



N 



N t'l 



(D17) 



(D18) 



(D19) 



where jiQand piare unknown mass scales determined by renormalization conditions. We set piQ = to have ( (c^)^) 
for flat space-time. 



3. Adiabatic Expansion of Time Integral 



In order to compute some correlators using the rn-rn formalism 1 32 1, such as two-pornt function {(t^C)' we need 
to integrate PV field contributions over time. In this subsection, we present how to calculate the time integral of PV 
fields by adiabatically expanding the integral. 

For simplicity, consider a diagram with one internal vertex. Using the WKB solution | D14| of a PV field, the 
general form of the time integral is 

I(/ci,/c2,- ■ ■ ,tf) = r dtGik^Mr- ■ ■,tf,t)e-'l'>^^^''^'*', 



(D20) 

where a;(f ) = o^j-j {t) + o^j-^ (f ) + ■ ■ ■ and G{ki, k2,--- 'AfA) = 0{co"). Because the integrand is a rapidly oscillatory 
function, the dominant contribution comes near the final time tjr. Thus, using integration by parts we expand the 
integral with respect to co: 



I{ki,k2,- ■■ ,tf) 



G{ki,k2, - ■ ■ -Jfjf) 



ico{tf) 
G{ki,k2,- ■ ■ ■,tf,tf) 



dt 



dt icv{t) 
1 d G{ki,k2,--- -JfJ) 



(D21) 



ILO 



1 



I ico{t) dt ico{t) 
1 d G{ki,k2,- ■ ■ ■,tf,ty 



t=tf 



iLo{t) dt \ ico{t) dt 



iLo{t) 



+ 0(a;"-*). 



(D22) 
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Note that the mode functions m„ j- and m* j. appear in pairs because of Wick contraction. Hence, the final result should 



be written in terms of 



Un,k{if) and their time derivatives, and we can compute the time integral up to arbitrary 



order of co. It is straightforward to generalize this to the cases with any number of internal vertices. 



Appendix E: Two-Point Function ( (cr^) ^ f ) in the Uniform Curvature Gauge 



In this section, we compute ( (y^),. V) using the uniform curvature gauge in the quasi-de Sitter(dS) background, 
where the slow-roll factor e is constant. Then we will show that the results in the both gauges are consistent with 
each other. Particularly, for the massless limit, the next leading order term in the uniform curvature gauge that 
indeed decays as / a^. 

The two-point function is the same as in the comoving gauge except that the counter term contribution appears in 
the leading order. 



(El) 







j d^xe- 






J N=0 \ 



where R is the Ricci scalar. After taking non-derivate interaction term T^^Sg^^^"^ only, factoring e and ^ out from the 
integral, we get 



e / d^za'it,) Y: 

N=0 



4(t,f),(T; 



+24eH^ 



<5Zi + O e, e^. 



.2 P 



where we have used the perturbed curvature in the uniform curvature gauge 

R = 12H^ - GeH^ + lAeH^^ + AeHt, + ■■■ 

where ■ ■ ■ denotes 0(e, e^) terms or terms proportional to the equation of motion of ^. 
Since TOO = £^ + ^ 



(E2) 



(E3) 



With TOO = £^ + 







N 


[(? 


1 


d 1 


Hdt \ 



f I" d4zfl3(t,)<J[4(f,x),4(z)]) 

, we have 



, together with the identities 193 ,196 1, and 



(E4) 



FN(f) 



N 



(E5) 



'N 



27r2 



|"N,/c| 



2M- 



'N 



1 d / o 



(E6) 



Although the rhs of Eq.dETJ is well-defined and regulator independent, individual terms are not. Thus, we insert 
counter terms to have each term regulator independent 



N 



e(^2( [a'{t 
1 
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where we have put the counter terms 5Zq and 5Z\R. into each one-point function, and the PV field contribution 



from the third term cancels with those from the other terms. Then, using the relation 1 75 1 one can find the rhs of 
Eq.( |E7| is consistent with the result | [99) in the comoving gauge in the quasi-dS backgroundafter explicitly computing 

renormalized one-point function [cp-if])^ On the other hand, the rhs does not depend on the renormalization 

as all counter terms cancel. Hence, we can arrive at the same conclusion using the one point function using super- 
horizon approximation in the dS space-time, 

^'"'"'i^X - 

where the arbitrary constant c < 0(1). Note that the UV boundary of the integral should be a comoving scale in 
order to to keep the spatial dilatation symmetry. 



Massless Limit 



For the massless limit rn^/H^ ^ \np/aH, we can compute the two-point function explicitly without neglecting 
any gravitational couplings. We calculate up to the next leading term here. We decompose Eq. ||E1| as 



n=l 



(E9) 



where Iq, In, and let. are the contributions from the physical field cr, the PV field Xi and the counter terms, respec- 
tively. Since all the gravitational couplings are 0{e) (See Eq. |C26 ), we may use the mode functions and mj- in the 
pure dS for 0{e) correction to the two-point function. Then a long but straightforward calculation gives 



J N=0 
1 p3 A 



+0{e^,e). 



TO 



A 5 



A 



3 fl3H3 flH 



+ 21og- + -4--2log 
p 3 a-'-H^ p 



O 



(ElO) 
(Ell) 



The PV field contribution I„ requires some more technical explanation. If we write the WKB solution IDII I as 
the PV field contribution I„ is written as 



where 



G„(fci,fc2;f,fz) = -2fl3fp(0afci(0a)c2W Eo?k;(fz)<(fz)4,(tz), 



Oi 



2 ^ ^^^n 



(fa;,^(fz)+a[P) (fa;,,(fz)+af ) 
^-P^^-...,(^.)+a«)+^(.a.,,(..)+3|f)) 



(2e), 



(E12) 
(E13) 

(E14) 
(E15) 
(E16) 



where d\' and 3^ denotes the time derivative with respect to a| {tz) and ^*p{tz), respectively, and \ 0\j and \ 0\j cor- 
respond to the (00) and the (iO) components of the gravitational couplings, respectively. Notice that aj- = O ^a;^^'''^^ 
and G (fci, k2; t,tz) = 0{lo^), and thus In has quadratic divergences superficially. However, the quadratic divergences 
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arising from (Oi^ vanish in the M„ oo limit. Effectively, the integral |E13 1 is linearly divergent. That means we 

have to adiabatically expand the integral to the second order. Similarly, the mtegral of the two-point function in the 
comoving gauge is quadratic divergent, and thus one need to expand the integral to the third order. This makes the 
computation easier in the uniform curvature gauge. Using 



1 



1 
1 



'yo{k,t)-icv^h{^,i)+0{cvf) 



h2 / „2 

-3iH - 2ijo{k, t) + i^^H + 0(a;-i) 



where 



3 IfcVa^ 
2 2 col 



H, 



H + 2H2 ^ (H + 3H2) Ml 5H^Mt 



which are obtained by combining Eq. | D14| with Eq. |E12 1, the integral | E13| becomes 



InipJ) 



TO 

bp 



1 A ^, 2A 5 p2 , 2A 



3 fl3H3 aH 



25 p2 



+ 



aM„ 3fl2H2 ^ aMn 3 
+ 0(ele). 



Note that all A depe ndent terms in Iq + Yj„ I„ vanishes by the PV field normalization conditions | |D3| 
Putting Eqs. | |E11| and |E22 1 together into Eq. ||E9|, we have 



TO 

bp 



am 5 p2 flf/i 8 
'l°§2^ + 3^1°§2^ + 3 



7 



18 fl2H2 



+ 0(e^e). 



(E17) 
{E18) 

(E19) 

(E20) 
(E21) 



(E22) 



(E23) 



We still need to compute one-point fimction ^ uptoO(e) in order to compare the results in both gauges. Be- 

cause mode functions for a massless scalar field are 0{e^), we need 0(e) correction on it. In a quasi-dS background, 
we take an ansatz for the mode function 



Uk{i) 



V2ka{t) 



+ i 



m 



ga{l)H{t) 



et) 



lka{t) 



(E24) 



where fk{t) = 0{e^) so that it recovers the dS solution in the e ^ limit. Applying this to the differential equation 



k^ 

%(f)+3HMfc(f) + -2"/c(0 =0, 



{E25) 



we get 



whose solution is 



a{t) 



fc2 



/, + H{t) - li—r h - H{tyfk = - 2/-77tH(0 - 2—. + 0{e 



3 . il 

-2+^^+2^ 



+ci 1 + 



C2 1 



a{ty 

e-^"lEi{2iq) 

2iq 



(E26) 

(E27) 
(E28) 
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where q = ^^jj^jyIJ)' exponential integral function 



Ei{z) 
Ei(±ix oo) 



dt 



±iTZ + e 



±ix 



0! 



1! 



2! 



(±fx) (±fx)2 (±fx)3 



(E29) 
{E30) 



Matching this solution with the Bunch-Davies initial condition | |D8| and the Wronskian condition | |D9| respectively 
give 



C2 = —in and ci = -. 
Then the mode function with 0{e) correction in a quasi-dS space-time becomes 



1 . H \ 

+ 1 ^ 



2ka 
e 

2ka 



/2Fy 

, . k .aH 



1-i 



.aH 



-in + Ei{2i^) ] e~^'M 



e'lH + 0{e^,£). 



Now we calculate the one-point function using this mode function as shown in Subsection D 2 and we get 

H 



dt 



H3 eH3 
4^ ^ 2^ 



log 



7 +0(e^e). 



Finally, we find 

1 d 
Hdi 



47t2 
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\ aH 3 
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»i4 



^-27+3log^+21og^ 



+0 e^^e. 



The non-p2 / fl2-suppressed terms are rewritten as 



47r2 





2 eH2 
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+ 27r2 
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log 2^ + 2-7 



47r2 
47r2 
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1 + 2e log 



TO 
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fJLV 
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(E31) 

(E32) 
(E33) 

(E34) 



(E35) 

(E36) 
(E37) 
(E38) 



As expected, this is the result (99 1 in the comoving gauge. The other terms are suppressed by the factor p / a . This 
explicitly proves that the next leading terms for the two-point function ( ((j2)^ are 0{p^ / a^). 
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